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We develop a Green's function based mean-field theory for coherent mixing of matter- and light- 
waves. To demonstrate the utility of this approach, we analyse a co-propagating Raman matter-wave 
amplifier. We find that for a given laser intensity, a significantly faster amplification process can be 
achieved employing resonant rather than off-resonance driving. The ratio of the matter-wave gain 
to atom loss-rate due to spontaneous emission is given by the optical depth of the sample, and is 
the same both on- and off-resonance. Furthermore, we show that for short-times, the single-mode 
approximation for the matter-waves gives exact agreement with the full spatial dynamics. For long 
times, the off-resonant case shows suppressed amplification due to a spatially inhomogenous AC 
Stark shift associated with laser depletion. This suppression is absent on-resonance, where the AC 
Stark shift is absent. 
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I. INTRODUCTION 

Matter-wave amplification (MWA) in atomic ensem- 
bles is a stimulated (i.e. collectively enhanced) nonlin- 
ear wave-mixing process in which two coherent matter 
waves (pump 1 + signal) and one coherent light wave 
(pump 2) are input, resulting in amplification of the sig- 
nal matter wave, along with the phase-coherent gener- 
ation of an additional light wave (idler), whose wave- 
vector is determined by phase-matching conditions [l]- 
[50]. As the amplification process preserves the relative 
phase between the two matter waves, MWA has been 
viewed as an active device with potential applications in 
matter- wave interferometry |21H26) . This point of view, 
however, is somewhat limited by the fact that the signal 
is amplified at the expense of pump 1, whose depletion 
eventually suppresses the MWA gain mechanism. On the 
other hand, when viewed as a fully quantum-mechanical 
process, a MWA-tyle process can in principle map an 
arbitrary quantum state from a two-mode atomic field 
(pump 1 -|- signal) onto a two-mode optical field (pump 
2 -I- idler) [3l HI [15] . Thus quantum information acquired 
by a matter- wave field could, e.g., be transferred to an 
optical field for detection or long-distance transfer to a 
second bimodal atomic system [27l [28] . 

As a mechanism for information transfer, MWA is com- 
plementary to electromagnetically-induced-transparancy 
(EIT) based quantum information retrieval techniques 
[SO]- Whereas EIT approaches use an adiabatic 
'counter-intuitive' pulse sequence to deplete the small sig- 
nal mode, while mapping the stored information onto an 
outgoing light pulse; the MWA process uses a dynami- 
cal 'intuitive' pulse sequence to amplify the small signal 
state, mapping a complementary quantum state, corre- 
sponding to the amplified signal state minus the origi- 
nal signal state, onto the outgoing light pulse. For the 
idealized case of 100% pump-to-signal conversion, this 
complementary state corresponds to the quantum state 
of pump 1. If the initial populations of the pump and 
signal matter-waves are equal, then the EIT and MWA 



information retrieval techniques would be effectively in- 
distinguishable. 

In developing quantum information protocols based 
on atomic systems, a universal challenge is to overcome 
the inevitable collisional dephasing and/or environment- 
induced decoherence processes. For MWA-based devices, 
one possible solution is to use very strong pumping fields 
to complete the process on a time-scale where dephas- 
ing/decoherence is negligible. In order to determine if 
this approach is feasible, we must first determine the ulti- 
mate rate limit to MWA. Most MWA and closely-related 
'matter- wave superradiance' experiments [H[5l [TTl[T?[[5T] 
operate in the regime of large detuning, where MWA can 
be attributed to collective Rayleigh or Raman scattering. 
It is clear, however, that the scattering rate increases as 
the detuning is decreased [TTl [50] , reaching its maximum 
value at zero detuning. With this in mind, we propose 
performing MWA on resonance, employing a short, in- 
tense pulse for pump 2, thus maximizing the underlying 
photon emission rate. 

Typically, the atoms used in MWA are cooled well be- 
low the critical temperature, Tc, for Bose Einstein con- 
densation, where the spatial coherence length reaches its 
maximum possible value, i.e. the sample size. For this 
reason, a BEG is an ideal source of coherent matter waves 
in much the same way that a laser is an ideal source of 
coherent light-waves. Despite these considerations, co- 
herent nonlinear mixing of light- and matter-waves does 
not require a BEG [3H [33] , and has been demonstrated 
experimentally for T > Tc [33]. In this work, however, 
we will focus solely on the BEG regime, with the under- 
standing that many of our conclusions may apply also to 
thermal g the two systems are indistinguishable 

for timescales small compared to the coherence lifetime, 
^ds/vr, where Xds is the thermal de Broglie wavelength, 
and Vr is the recoil velocity. 

While going on-resonance is often associated with 
rapid heating, and matter-wave decoherence, this is not 
actually the case in MWA. While going on-resonance cer- 
tainly increases the decoherence rate, it also increases the 
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MWA gain rate. As both of these rates scale linearly 
with the background emission/scattering rate, their ra- 
tio, which determines the MWA fidelity, is a universal 
property, i.e. independent of the detuning and pump 
strengths for a wide range of parameter settings. This 
means that a short, intense control pulse (pump 2), could 
complete the MWA process on a time-scale much shorter 
than the spontaneous lifetime, with the net heating ef- 
fects being no worse than those found in off-resonance 
MWA. 

This approach cannot be applied to Rayleigh scatter- 
ing MWA schemes, where reabsorption of the idler pho- 
tons will lead to a cascade of many hundreds of ad- 
ditional matter- wave momentum sidemodes llj on the 
MWA timescale. It therefore applies only to Raman 
MWA, where a different set of problems arise. In the 
case of Rayleigh scattering, the AC stark shifts imposed 
by pump 2 are the same for the signal and pump 1, so 
that the relative phase between the two matter waves re- 
mains constant. In the case of a Raman scheme, however, 
the AC stark shifts imposed by pump 2 do not match, as 
different hyperfine states generally have different dipole 
moments. As a result, the matter-wave relative-phase 
will oscillate rapidly during the MWA process (typically, 
the oscillation rate will be ^ 100 times faster than the 
MWA gain rate). This effect, however, does not occur 
in the resonant case, where the AC Stark shift vanishes. 
Thus in addition to providing the largest possible gain, 
resonant MWA also eliminates the undesirable AC Stark 
imbalance inherent in Raman MWA 

Ultimately, the source of dynamical nonlinearity in 
MWA is the atom-field dipole interaction, which becomes 
a cubic non-linearity upon second-quantization of the 
matter-wave field. The interaction between the matter 
and light fields is primarily governed by four important 
parameters: the Rabi frequency, SI, the detuning A, the 
spontaneous emission rate, F, and the optical depth, D, 
which governs the strength of the nonlinearity. 

In the far-detuned (|A| ^ |r2|,r_D), and overdamped 
(A = 0, TD ^ \n) regimes, the excited atomic state 
is only virtually populated, and can be adiabatically 
eliminated. This results in an effective cubic nonlin- 
earity whereby MWA becomes a four-wave mixing pro- 
cess, unique only in that the waves are of two different 
types (2 matter -I- 2 light). In the Rabi regime (A = 0, 
|S1| S> ri?), however, the excited state is macroscopically 
populated, so that the system consists of five dynam- 
ically independent waves, which interact via a pair of 
coupled cubic nonlinearities. Thus, aside from possible 
practical applications, resonant MWA may be a system 
of fundamental interest from a nonlinear-optics perspec- 
tive. In fact, treating the excited level dynamically is 
equivalent to extending the nonlinear susceptibility to all 
orders, with the Rabi regime corresponding to the non- 
perturbative limit of a saturated absorber. In practice, 
the resulting strong nonlinearity can lead to nontrivial 
spatial effects in both light and matter waves, as was 
observed in a recent Rayleigh superradiance experiment 



under near-resonant driving 



A. Preliminary estimates 

Emission of a photon imparts a recoil-kick onto the 
emitting atom. The solid angle of emission such that 
the final momentum of the atom lies within the atomic 
momentum coherence length of a given direction defines 
a Fresnel mode. Even in a thermal sample, the atomic 
momentum coherence-length is non-zero due to the finite 
size of the sample. For a thermal atomic vapor of width 
W, the corresponding thermal momentum coherence- 
length is kcoh ~ ^/W, so that the Frcsnel-mode solid 
angle is iljr = TT{kcoh/kLf n/ikLW)"^. The MWA 
idler mode corresponds to a Fresnel mode aligned in the 
direction for which an atom excited from pump 1 is trans- 
ferred by photon-recoil into the signal mode. In the ab- 
sence of signal population, the probability of emission 
into this mode is approximately flj^/An <^ 1 (neglect- 
ing dipole emission pattern effects). Bosonic stimulation 
should enhance emission into this mode by the signal- 
mode population, {Ns + 1), relative to all other emission 
channels, so that the probability for emission into the 
idler mode is 



P,. = 



4n + Nsfljr 



(1) 



with non-idler emission occurring with probability I — Pi. 
The ratio of idler to non-idler emission, which corre- 
sponds to MWA gain-to-loss ratio, is thus given by the 
collective-emission optical depth, 



D = {Ns + l 



47r 



(2) 



where Aj_ is the cross-section of the atomic sample per- 
pendicular to the idler direction. Thus a necessary condi- 
tion for MWA is a large collective-emission optical depth, 
Ds ^ 1- In terms of the density n = ns/{LW'^) this gives 
the standard result D ~ nX^L. Clearly an elongated, 
"cigar-shaped" sample of length L and width W, will 
maximize D for fixed atomic density, provided the signal 
mode is directed along the long-axis. Thus, at present 
we consider only condensates with aspect ratios greater 
than or equal to unity, i.e. L > W. 

Any Fresnel mode within the geometric solid angle 
Hef = tt{W/L)'^ of the long-axis is termed an 'end-fire 
mode'. Any mode outside this solid angle has an effective 
optical depth reduced by the factor W/L relative to that 
of an end-fire mode. The number of end-fire modes is 
then given by M ~ Hef/^j^ ~ klW^/L"^. This leads us 
to an important parameter, the 'Fresnel number', given 
by = k^W^/L, so that M ^ J^^. A typical cigar- 
shaped BEC, with an aspect ratio of 10 : 1, has ~ 800, 
which does not necessarily correspond to the optimal op- 
tical depth for collective emission. For fixed atom num- 
ber, the optical depth would be increased by increasing 
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the aspect ratio. In doing this, however, one eventually 
reaches the regime where J- < \. At this point, light 
emitted into an end-fire mode diffracts out of the con- 
densate volume before reaching the far end of the BEC, 
resulting in a decreased effective optical depth. This ef- 
fect is characterized by the diffraction length, Ld = k^w"^ , 
which is the distance a coUimated beam of light will prop- 
agate before it begins to spread. 

We see that the ratio of the diffraction-length, Ld to 
the sample length L leads to the Fresnel number T = 
k^W^/L. Thus for 1 <^ J^, the optical depth is sample- 
length limited, D ^ nX^L, which can be re-expressed 
as 



N r- 
lim D =V^, 



(3) 



where v = k^LW^ is the dimensionless sample volume. 
For ^ 1, the optical depth is instead diffraction-length 
limited, D ^ nX^Ld, which can be re-expressed as 



N 

Urn D ^ 



1 



(4) 



Thus we see that the maximum MWA gain for fixed atom 
number N is obtained by decreasing the volume while 
adjusting the geometry such that ^ 1, resulting in a 
maximum optical depth of 



N 

7^' 



(5) 



The per-atom rate of non-collective photon emission 
from an atomic sample is given by PeT, where Peit) is 
the single-atom excited state probability. Generalizing 
this to the nonlinear regime of collective emission into 
the idler mode, the resulting exponential gain rate for 
the signal mode occupation is then 



P.TD. 



(6) 



where Pe is the single-atom excited state probability, and 
r is excited-state spontaneous emission rate into the elec- 
tronic state of the signal matter-wave. To good approxi- 
mation, the excited state probability of a driven atom in 
this system is 



1^21' 



4A2 + {rD)^ + 2\n\ 



(7) 



where A, and Q, are the detuning and Rabi frequency, 
respectively, of pump 2. 

By inserting ([t]) into ^ , we can see that for fixed laser 
intensity, the MWA gain is always maximized by going 
on resonance, A = 0. Equation (3) shows that the atomic 
transition saturates in the limit jfip ^ ^{TDs)^, leading 
to the condition 



Gs < TDs/2, 



(8) 



as an initial estimate for the upper-limit to the matter- 
wave amplification rate. It is important to note, that the 



saturation intensity has increased by D^, the square of 
the collective-emission optical depth, due to superradiant 
broadening of the excited-state linewidth. 

The idea to use forward scattering of short, intense, 
resonant pulses to probe a BEC was first proposed in 1995 
by You, Lewenstein, and Cooper, [33], who considered 
only two-level atoms without recoil in the single spatial- 
mode approximation. While this paper established the 
feasibility of resonantly probing a condensate on a time- 
scale short compared to the spontaneous lifetime, show- 
ing that meaningful phase information could be obtained 
without decoherence, to the best of our knowledge, no 
experiments have been performed along these lines. Our 
proposal for resonant MWA is, from a formal viewpoint, 
an extension of this idea to the case of a lambda-type 
three-level system. 



B. System-reservoir versus Maxwell-Schrodinger 
models 



There have been two basic models used to describe 
the MWA system. The first, pioneered by Moore and 
Meystre [6] , is based on extending the Wigner-Weisskopf 
theory of spontaneous emission to the regime of atomic 
field theory. In this approach, the light field is elimi- 
nated, resulting in a set of quantum Langevin equations 
for the atomic field [36]. These equations are further 
simplified by quantizing the atomic field onto a quasi- 
complete set of 'momentum side-modes' [51 [371 [35] ■ We 
will henceforth refer to this approach as the quantum 
coupled mode (QCM) model. In this model, the MWA is 
attributed to bosonic stimulation by the atoms, whereby 
recoil into a particular side-mode is enhanced by a factor 
of the side-mode occupation number. 

The QCM mode, however, fails to explain several im- 
portant observations [TTl [SH]- These include asymme- 
try between forward and backward superradiance modes, 
spatially dependent depletion of the BEC, and subexpo- 
nential growth and ringing of the scattered field [101 131] • 
These effects, usually referred to as spatial effects or 
propagation effects, are shown to be important in the 
strong driving regime, even though atomic motions can 
be neglected [H]. In this case, whereas the total atomic 
density is a constant of motion, the distribution of prob- 
ability over internal energy levels can still depend on po- 
sition. 

In a more recently developed approach aimed at ad- 
dressing these phenomena, MWA is viewed as Bragg 
scattering process with feedback, in which atoms are 
diffracted by optical standing waves formed by the pump 
and idler fields [20l l40ti47] . The superradiant enhance- 
ment occurs because the diffraction efficiency is propor- 
tional to the square to grating depth, which depends 
on the atom number in the signal wave. This ap- 
proach incorporates the spatial dependence of the light- 
matter coupling via coupled Maxwell-Schrodinger equa- 
tions, which can be efficiently solved at the mean-field 
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level, hence we term this the semiclassical Maxwell- 
Schrodinger (SMS) model. As the initial build up of the 
superradiance from noise is not a mean- field effect, it was 
treated in the QCM model, which was used to generate 
suitable initial conditions for the SMS equations. For the 
case of MWA, however, there is no build-up from noise, so 
that the dynamics is essentially coherent, thus mean- field 
theory is expected to work extremely well, as phase fluc- 
tuations are effectively negligible for large enough parti- 
cle numbers. 

At the mean-field level, and in the single-spatial mode 
approximation, the QCM equations reduce to three or 
four dynamical variables, so that the dynamics can be 
well understood, and accurately described by approxi- 
mate or exact analytic solutions. Thus we will first ex- 
plore the QCM results, and then compare and contrast 
them with the more accurate, but numerically intensive 
SMS dynamics. 

We emphasize that in both models, the light fields adi- 
abatically follow the matter-wave fields, the difference is 
that only the SMS respects causality, in the sense that 
the intensity at position r of a field mode propagating 
in direction k can only depend on the integrated atomic 
polarization in the —k direction. The point is that one 
should not be misled into thinking that the light fields are 
independent dynamical variables in the SMS approach. 
In general, the coupled atom-field equations can always 
be replaced by non-local, but purely atomic non-linear 
equations using Green's function methods. 



C. Organization 

In order to compare the on-resonance dynamics with 
the previously studied far-detuned dynamics, we will 
study the MWA process in the strong-pumping regime 
where pump laser depletion is a small effect. We will 
therefore explore the cross-over between two limiting 
cases: the Stark regime, |A| ^ |ri|,rD, and the Rabi 
regime, ^ FD, A. In order to illustrate important 
underlying physics, we will focus on a simple, three-level 
Raman system, with co-propagating pump and idler light 
fields. 

In section [ll A[ we use a modified master-equation ap- 
proach, treating the optical field as a Markovian reser- 
voir, and derive a set of non-linear non-local equations 
for the atomic field operators. Then in section II B[ we 
use this approach to derive the mean-field theory for 
both atomic and optical fields. Using a Green's func- 
tion formalism, we show that, at the mean-field level, 
treating the field as a Markovian reservoir is equiva- 
lent to solving coupled Maxwell and Schrodinger equa- 
tions, with the additional benefit of including vacuum- 
fluctuation effects such as spontaneous emission and the 
Lamb-shift. We will then compare a theory using the 
exact Green's function of an ideal dipole with previous 
approaches to describe MWA and superradiance in the 
Maxwell-Schrodinger picture. 



In section III we develop the quantum coupled-mode 



model, based on a single spatial mode approximation for 
the condensate, and give a precise calculation of the opti- 
cal depth and dipole-dipole shift . Then, in section [TV] we 
solve the full spatial equations in mean-field theory, in- 
corporating laser depletion and field-propagation effects, 
effects. Lastly, we give a discussion and conclusion in 
section |Vl 



II. THE MODEL 

A MWA process can be implemented with Raman or 
Rayleigh transitions, corresponding to the pump I and 
signal matter waves in different hyper-fine and/or trans- 
lational momentum states. In a Rayleigh MWA process, 
multiple atomic recoil modes, in addition to the initial 
pump and signal, would unavoidably appear due to two- 
photon Bragg scattering of atoms by the pump and idler 
fields [iTl |40l [42] . In a Raman MWA, however, such pro- 
cesses are limited by the number of ground-state hyper- 
fine Zeeman sub- levels [371 SI] ■ For this reason, Raman 
transitions are better suited to MWA applications, such 
as matter-to-light quantum state transfer. 

For concreteness, we consider a very simple Raman 
MWA configuration, with the pump laser and idler field 
corresponding to co-propagating right- and left-circularly 
polarized fields, respectively, i.e. a co-propagating A- 
system. In this configuration both the pump and sig- 
nal condensates are at rest (k fa 0). Thus the two con- 
densates must correspond to two different hyperfine mp 
states, separated by Amp — 2. This co-propagation ar- 
rangement is different from the commonly used end-fire 
configuration where pump 2 and idler beams are per- 
pendicular O IllL 115) . Comparatively, the pump-laser 
depletion effect is thus expected to be stronger in our 
scheme. 

As we will see in Sec. 



Ill B the maximum obtainable 



optical depth is given by D ^ 3.3 X N/y/kfy, where 
N is the atom number, and V is the volume. For a 
typical condensate with N ^ 10^ and k^V ^ 10^, this 
gives D ^ 10^. The superradiance-broadened lifetime of 
the excited state is given by TD. For a typical atomic 
linewidth of F ~ lO'^Hz, this results in TD ~ 10^°, which 
is larger than the typical excited state hyperfine splitting, 
Ahf ^ 10^. For an alkali metal condensate, which has a 
single valence electron with a ^6*1/2 ground-state config- 
uration, absorption of a (t+ photon and subsequent emis- 
sion of a a~ photon results in a hyperfine angular momen- 
tum change of Amp = 2. As hipping the electron spin 
can only give Amp = 1, the co-propagating A-system re- 
quires a change in the nuclear spin of at least Amj = 1. 
As the optical fields do not efficiently couple to spins, this 
nuclear spin-flip can only be accomplished by the hyper- 
fine interaction when the atom is in the electronically 
excited state. This therefore requires TD <C Ahf, i.e. 
D <C 10^. Under such a constraint, incoherent radiation 
and the corresponding loss of atoms from the pump and 
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FIG. 1: (Color online) Raman MWA system configured on 
a metastable J—1 state for an isotope with no iiyperfine- 
structure. Figure (a) shows the A-level scheme, with the laser 
and idler fields indicated. The pump condensate is in internal 
state |1) and the matter- wave signal is in internal state |2). 
The pump laser and idler fields are indicated by their respec- 
tive Rabi frequencies, Ql and fi^. Figure (b) illustrates the 
schematic setup for copropagating beams. 



signal modes will no longer be negligible. Therefore the 
standard alkali metal atom condensates are not suitable 
for resonant MWA experiments. 

The solution is to choose an atomic species with no nu- 
clear magnetic moment, / = 0, and therefore no hyper- 
fine structure. To obtain a A-system with co-propagating 
circularly polarized fields thus requires a ground state 
with J > 1. To avoid sequential re-scattering of the 
pump and idler fields, an ideal configuration would be 
J = 1, with the pump and signal matter- waves corre- 
sponding to rrij — —1 and mj = I. The level diagram 
and setup for a MWA scheme based on these considera- 
tions is shown in Figures [l^ and[T]D, respectively. 

To date, the only atomic species that has been con- 
densed into a J = 1 state is "^He* [IHl EH] , which has a 
■^Si configuration. Unfortunately, this choice is not vi- 
able because the excited-state fine-structure splitting is 
anomalously small at Afg ~ lO^^Hz. The best possi- 
bility is therefore to prepare a BEC with a species that 
can be coherently driven into a metastable J—1 state. 
Two atomic species with suitable metastable states for 
which BEC has been demonstrated are *°Ca [50] and ^''Sr 
[511 152] , both of which have zero nuclear magnetic mo- 
ment (/ = 0). The 4s4p level of ^^Ca has a hfetime 
of ~ Ims, and the 5s5p '^P° level of ^"^Sr has a lifetime 
of ^ 10/is, both of which would be adequate for resonant 
MWA, which operates on the time-scale of ~ Ins. An- 
other promising species would be ^'^^Ba, which has been 
laser cooled and trapped [53], with good prospects for 
achieving BEC [54]. Decay of its 6s5d level, which has 
a ^Di 2,3 configuration, is doubly-forbidden, and has a 
measured lifetime of ~ lOs. For all three of these species. 
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4s4p 3p° 



4s2 ISn 




6s5d 3D 



6s2 >Sn 



FIG. 2: (Color online) Two viable implementations, one using 
a metastable triplet P state in *°Ca and the other using a 
metastable triplet D state in ^''^Ba. By driving a J = 1 to 
J = 1 transition, a pure A system is realized. 



the fine structure splittings are ~ 10^'^Hz, which makes 
them ideal for resonant MWA experiments. In figure [2] 
we show level structures for potential MWA schemes con- 
figured on metastable J—1 states in ""^Ca and ^■^^Ba. 
By driving a transition to a J = I excited state, there is 
no coupling of the pump laser to state |2) nor of the idler 
field to state 1 1) , making the system a pure A configura- 
tion. 

While a Bose-Einstein condensate (BEC) is not re- 
quired for coherent collective effects such as MWA or 
superradiance, we will assume a BEC as a source due to 
its high optical depth. Later, we will consider the pre- 
cise temperature requirements imposed by each dynami- 
cal regime. In the present scheme, both pump and signal 
modes are at rest, while the excited state propagates at 
the recoil velocity. We presume this will improve the 
MWA fidelity, which will be important for applications 
such as state mapping between matter and light waves. 
We note that in our dynamical simulations, we will ne- 
glect the kinetic energy operator for all states, as well 
as s-wave collisions, based on the assumption that these 
only become important on time-scales much longer than 
the fast MWA time-scales we are interested in at present. 
The effects of s-wave collisions are included only via the 
initial Thomas-Fermi wavefunction of the BEC. 



A. Quantum-field theory for laser driven atoms 

The mean-field equations for the atomic field in the 
QCM model can be derived by treating the electromag- 
netic vacuum as a Markovian reservoir. For a generic 
system, governed by the Hamiltonian, Hs, coupled to a 
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zero-temperature bosonic reservoir, governed by 

Hr dkWkflkA^kA, (9) 

with a coupling described by 

y = ^ / dk (4,akA + aL^kA), (10) 

A 

where {skA} is any set of system operators. 

We proceed by separating the system operators, SkA 
into fast and slow parts, relative to the reservoir relax- 
ation time-scale, as 



SkA 



E 



(11) 



where SkAm is slowly-varying, and all fast dynamics is 
accounted for by the exponential factors. Note that we 
are not making the rotating-wave approximation, which 
would involve keeping only terms for which = 0. 

The equation of motion for the expectation value of an 
arbitrary system operator, S, is then given in the Markoff 
approximation by 

dt{S)=i{[H,+Vs,S])+Y, J dkirSiio^ - 

Xni 

^ (2SkAm5'SkAm — sj^.^^SkAm'S' — S'Sk^j^SkAm) , (12) 



where 



dk 



Wk — w„ 



■^kAm^kAm, 



(13) 



is the system potential induced by the interaction with 
the reservoir. 

The Hamiltonian for a system of three-level bosonic 
atoms coupled to the EM-vacuum, and driven on the 
|1) — > |3) transition by a laser field, Ei(r, t) — 
Ei(r)e~*"^* + E2(r)e*"^*, is given in the rotating frame 

by 

Hs = A ydr*J(r)*3(r) 

dr [d* • EL(r)*l(r)«'3(r) + i/.c] , (14) 

where A — — ujl is the laser detuning, and dj is 
the dipole moment of the \j) — >■ |3) transition divided 
by h. The atomic field operators obey the commutation 
relation [^'^(r), «',v(r')] = %'5(r - r'). 

The reservoir is characterized in the rotating frame by 
mode frequencies ojk — c{k — kL), and corresponding sys- 
tem operators 



*kAl 



E 

3=1:2 



d* ■ ekAj 



ck 



2eo(2 



(15) 



and 



*kA2 ^ E 

3 = 1.2 



Si^^ = y (dj • GkA) 



(16) 

with corresponding oscillation frequencies lji — Q and 
L02 = — 2cfci. The unit- vectors {ekA}; A = 1,2 are any 
two orthonormal polarization vectors that are mutually 
orthogonal to k. 



B. Green's function approach to mean-field tlieory 



Using (12 1 to compute the equation of motion for 
i!j{v) — (5'j(r)), and then making the usual mean-field 
factorization, we obtain the mean-field theory: 



idi ■ E* ^3, 

id2 • E* tAs, 
■ 7' 



-03 



-Ed.* 

3=1,2 



(17) 
(18) 

E0„ (19) 



where 7 is the excited state spontaneous decay rate, and 
E(r) is the positive-frequency part of the optical field. 
The spontaneous decay term in Eq. ( 19 ) accounts for the 



incoherent radiation emitted by the atoms, and results in 
an irreversible loss of the mean-field density. The optical 
field, E, is then the sum of the laser field and the coherent 
radiation emitted by the atoms. 

The total optical field is the sum of the laser field and 
a scattered field, E = E^ + Es, where the scattered field 
is related to the positive-frequency polarization density, 

P(r) = 5]d,V;(r)V3(r), (20) 



via a Green's function tensor, G(r), according to 

EJr) = ^ / dr' G(r - r') • P(r'). (21) 
eo J 

Let Go(r) be the Green's function tensor inferred from 
the equations of motion for the atomic field operators 
obtained via Eq. (IT2| . We find 



Go(r) 



2(27r)3 



dk(l-kk) 



2e_ 



- inS{k—l) 



(22) 



where ^ = k^r. Here kk is the (outer) tensor-product 
of the unit-vector k with itself. Performing the /c-space 
integration leads to 



Go(r) 



+ 



2 



1— f f + (1—3 f f ) 



1 

e 



(23) 



7 



The field generated by contact term in Eq. ( |23[ ) cor- 
responds to the internal field of the atom in the ideal- 
dipole limit |55| . To form a proper mean-field theory, we 
must drop this contact term, as an atom can feel neither 
its own internal field nor the field inside another atom. 
Subtracting the internal fields of the atoms is known as 
the 'local- field correction', and has been shown |551 [55] 
to lead directly to the Clausius-Mossotti/Lorentz-Lorenz 
relation. The proper Green's function is therefore given 

by 



G(r) = 



Mr) - ^'^'W, 



(24) 



which leads to 



G(r) = ^ /rfk(l-kk) 



fc2-l 



Performing the fc-space integration then gives 



1 — ff+(l — 3rf) 



i 1 



(25) 



(26) 



which reproduces via E(r) = 'G(r) • P, the well-known 
electric field of an oscillating ideal dipole located at the 
origin [571 . 

This particular formulation of the Markoff approxi- 
mation includes retardation effects related to the fast- 
oscillations of the atomic dipoles at the laser frequency, 
but neglects retardation related to the slow dynamics of 
the atomic fields in the rotating frame. This is valid as 
long as the time-scale of the slow-dynamics is short com- 
pared to the retardation time-scale tj. — L/c, where L is 
the system size, and c is the speed of light. For a typical 
BEG, we have L ^ 100 /zm, so that tr 1 ps. In contrast, 
the fasted time-scale in the rotating frame is given by 
the super-radiance-enhanced decay rate YD. For typi- 
cal optical depths of 13 ^ 10"^, this gives a time scale 
of r ^ 100 ps. For larger systems, and/or higher opti- 
cal depths, the theory can be reformulated to take into 
account all retardation effects. 

The mean-field equations of motion, ( ITpQ I, can be 
expressed as a set of non-local non-linear equations for 
the matter waves, given by 

5*^1 = + \ J2 Gl.m^A^s. (27) 



3 = 1.2 



(28) 



5*^3 = -*[A-Z^]^3 



l^L^^i + l E GjkmHi^j, (29) 



where flLir) = 2d* • EL(r), 
2fc? 



G.,(r) = ^d:.G(r).d„ 
Co 



(30) 



and we have introduced the compact notation 

G[/]= /dr'G(r-r')/(r'). (31) 



The re-scaled Green's functions, which now have units of 
s"-'^, are given For di = ^(x + iy) and d2 = ^(y + ix), 

by 



G«(r) 



3r 



,.2 ' \ "^^2 "'^ M ^ ^2 



Gi2(r) 



and 



G2i(r) 



.'3Tix-iyy A 3i 3\ e'« 



. 3r (a: -f iy) 



3i 



with 



r = 



3vreo 



(32) 
(33) 

(34) 

(35) 



The diagonal fields, Gn and G22 describe two-body in- 
teractions that conserve atomic hyperfine angular mo- 
mentum. The off-diagonal fields, G12 and G21, describe 
two-body interactions that convert atomic hyperfine an- 
gular momentum into center-of-mass orbital angular mo- 
mentum, with only the total angular momentum being 
conserved. Only the diagonal fields contain 1/r terms, 
and therefore dominate in extended systems. The initial 
matter wave signal and pump fields are configured so that 
the MWA process is driven by these diagonal fields. The 
off-diagonal fields then describe the generation of vortices 
during the MWA process. The effects of the off-diagonal 
fields, however, are expected to be extremely small com- 
pared on the MWA time-scale, so that these processes 
can be ignored. 

Treating the off-diagonal fields as a perturbation, and 
factoring out the rapid spatial oscillations of the excited 
state and the driving laser via '03 (r) = 0e(r)e*'^^^, and 
r2i(r,i) — il(r, t)e*''^^, the unperturbed system evolves 
according to 



-(r!* + G*[v:0i])Ve, 



9*01 

5*02 = ^G*[V'*02]0e, 
9*0e 



0e 



(36) 
(37) 



where 



G(r) 



2k^ 



i=i,2 



(38) 



i 1 



(39) 

This set of equations will form the basis for our theo- 
retical description of the MWA process for the present 
configuration. 
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C. Comparison with previous approaches 

Previous theoretical models [40l |41] have described 
spatial propagation effects using coupled Maxwell and 
Schrodinger equations. The advantage of the Markovian 
approach is that it readily incorporates effects due to 
electromagnetic vacuum fluctuations, e.g. spontaneous 
emission. Aside from spontaneous emission, we will now 
demonstrate that in principle, the two approached give 
identical results at the mean-field level, the approxima- 
tion used by previous authors misses fundamentally im- 
portant physics. 

In the Maxwell-Schrodinger approach, the matter fields 
of a two-level atom are governed by 

dtiPg{r,t) = zd.E*(r,t)V'e(r,<), (40) 
dtipe{r,t) = -iWegVe + «d* • E(r,i)V'3(r,i), (41) 

while the electric field is determined by 



E(r,t) 



d^^;(r,t)^e(r,i)- 



(42) 



In order to solve this equation, one can first make the 
slowly varying envelope approximation, via 



E(r,i) 

V'e(r,0 



E(r_L,z,t)e*'=^(^""*), 



(43) 
(44) 



Keeping only the leading order terms gives the standard 
Paraxial Wave Equation (PWE): 



a. 1 - vi 



E 



/r2 



(45) 



For small-enough systems, we can neglect retardation, in 
which case the PWE reduces to 



[-2zfci9, - Vi] E = d^V^V', 



(46) 



For a system which is very slowly varying in the perpen- 
dicular direction, we can drop the transverse Laplacian, 
in which case we obtain the Longitudinal Wave Equation 
(LWE), 



2ikLd^E 



(47) 



which yields the equal-time relation 



E(ri,z) 



+ i 



E(r_L,2:o) 

dfcL 



2eo 



dz'rn{r±,z')Mr±,z'), (48) 



which was used to determine the electric field in previous 

works gnun]. 

Both the PWE and LWE can be solved using the 
Green's function formalism. The Paraxial Green's func- 
tion tensor satisfies 



(49) 



and is given by 

Gp(r) 



u{z) 
Attz 



exp 



2 1 



kLri 



(50) 



where u{z) is the unit-step function. With the definition 



G(r) = ^d* • G(r) • d, this becomes 



iTu{z) 
Gp(r) = exp 



2z 



2 kLZ 

The longitudinal Green's function must satisfy 
-2zfci9,GL(r) =<53(j.)^ 



and is given by 



so that 



GL{r)=i^u{z)S\v^). 



(51) 



(52) 



(53) 



(54) 



We note that the longitudinal Green's function is purely 
imaginary. This means that it can describe collective en- 
hancement of spontaneous emission (superradiance) but 
it does not produce the nonlinear dipole-dipole shift (i.e. 
the effective two-body interaction potential) . The parax- 
ial Green's function, on the other hand, has both real 
and imaginary parts, and therefore captures the essen- 
tial physics of the light-matter interaction, differing from 
the full Green's function only in its degree of precision. 
The Paraxial Green's function can be derived as an 



approximation to the full Green's function (39) by ex- 



panding separately the pre-factor and exponent to lead- 
ing order in powers oi rj_/z. which gives 



G(r) 



3F 
2kL\z 



exp 



2 1 



kLrj 
' 2|z| 



JkL{\z\-z) 



(55) 



Note that the phase factor e*'^^'^'^!"^) becomes unity for 
z > and e"*^*^^^ for z < 0. Thus, under integration we 
have 



JkLi\z\-z) 



u(z). 



(56) 



so that we obtain the paraxial Green's function of Eq. 
(51) as the leading-order expansion of G(r) is powers of 
r±/z. Thus we should expect the paraxial Green's func- 
tion to give good accuracy for condensates with large 
aspect ratios, and hence small J^, with the accuracy de- 
creasing as increases. 

To obtain the Longitudinal Green's function, we as- 
sume that the polarization density is approximately uni- 
form on the scale of the optical wavelength, and replace 
the Gp(r) and replace the transverse exponential with a 
transverse delta function, 



exp 



2z 



(57) 
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where 



j drj^ exp 


I — 




2z 



2ttz 



(58) 



which leads directly to Eq. ( 54 1 



In section |IV[ we will compare the mean-field results 
for the full Green's function to those given by the Parax- 
ial and Longitudinal approximate Green's functions and 
the single spatial-mode approximation presented in sec- 
tion Hill 



III. QUANTUM COUPLED-MODE THEORY 

In this section, we develop the quantum coupled-mode 
(QCM) theory for the three-level Raman MWA system. 
In section III A, we take the mean- field equations ( 36p8 ) 
and expand them onto a truncated set of spatial modes 
defined by the Thomas-Fermi (TF) density profile of 
a condensate. In section III B[ we then evaluate the 



complex collectivity-parameter which governs the MWA 
process. Because the TF density is a very accurate 
description, we expect quantitative agreement between 
the calculated parameter and experiment. In section 
[III C[ we develop approximate analytic solutions to the 
QCM equations of motion, for each of the three main 
dynamical regimes, which are then compared with exact 
numerical solutions of the QCM equatons in section 

HlKTsl 



A. Spatial mode Expansion 



In practice, the mean- field equations (36|38) are diffi- 
cult to solve numerically, due to their non-local nature, 
and the singularities in G(r). For short times, it may 
be valid to neglect the spatial dependences of the mean- 
field parameters and invoke the single-mode approxima- 
tion. This is accomplished by introducing the mode am- 
plitudes, Ci, Cg, and Cg, via 



V'i(r,t) = 
V'2(r,t) - 
V'3(r,t) = 



Nci{t)(t>{v) 



(59) 
(60) 
(61) 



where N is the initial atom number, 4'{t) is the ini- 
tial wavefunction of the condensate, and k/^ = ki^i the 
wavevector of pump 2. The atomic population in the 
pump mode is then Ni{t) = N\ci{t)\'^ , and the signal 
mode population is Ns{t) = N\cs{t)\^ . 

After inserting the mode expansion, the mean-field the- 
ory reduces to 



dtci = '-nic, + h:rN\c,\^c^ 
dtCs = ^r/*iv|ce|'c. 



(62) 
(63) 



dtCe = -i(A - i"f/2)Ce + - flLCl 



Cj I Cp. 



(64) 



where the collectivity parameter, /, is given by 

— / d^rd\'U(r)\^U(r')\^G(r -r')e-'''^'-^~^'\ 

«r J 

. 

We note that have has both a real and an imaginary 
part. The real part, fa leads to enhanced decay of the ex- 
cited state, with commensurate feeding of the two ground 
states. The optical depth is therefore given by D = JrN. 
The imaginary part, // leads to cross-phase modulation 
between the excited and ground states, thus imprinting a 
nonlinear phase-shift onto the signal mode. Introducing 
X as the ratio between the amplitudes of the real and 
imaginary parts, we have, 



f = fB{l-ix)- 



(66) 



With the dimensionless variables r = TDt, Q, — fi/FD, 
7 = j/^D, and A — A/TD, the equations of motion 
become 



drCi 
drCs 
drCe 



^^*Ce + ^{l+ix)\Ce\^Ci 

^(l + ^x)|cepc, 
-i{A - i7/2)ce + ^nci 



(67) 
(68) 

(69) 



These equations form the basis of the quantum coupled- 
mode (QCM) model. 



B. Collectivity parameter 

Within the QCM framework, the strength of the non- 
linearity is governed by the collectivity parameter, /, 
which depends only on the condensate geometry. By 
combining Eqs. (25) and (65), we can express the col- 



lectivity parameter as 
3 



f = 
where 



dk(2-fcl) 



5(k) 



7r(5(fc2-l)- 



fc2-l 



dV|0(r)pe 



S{kLik-z)), 
(70) 

(71) 



is the static structure function of the condensate. 

The spatial profile of the condensate is generally given 
to excellent approximation by the Thomas-Fermi distri- 
bution. 



15 



rW^L 



1 - 



(72) 
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where W and L are width and length of the conden- 
sate. From Eqs. (72 1 and (71 1, it is straightforward 



behavior, 



to show that for the Thomas-Fermi density profile, the 
static structure function is given by S'(k) = |p (iir(k))|^, 
where 



p{K) = 



15 



45 



45 

sin(i^) + — cos(if), (73) 



and 



K 



kyW 



k,L 
2 



(74) 



It is logical to parametrize the condensate dimensions, 
L and W, by the Fresnel number, T = ki^W^/L, and 
the dimensionless volume = k\LW^ , so that we can 
consider changing the geometry of the TF density profile 
while holding the mean-density fixed. With (k^W)'^ — 
\J vT and kj^L ~ -^/w/J-", and with the change of variables 
from {fcz,fc_L} — > related by k±^ — ^^^ x\/l~v? ^ 

and fcz — 1 = 2\ —xu, this becomes 



/ 



dx \p{x)\^ {^g{v,F,x) — ih{v,T,x)), (75) 



where 



and 



TrV^-2x^{l-ul) 
y (u+-u_) 



(76) 



h(v, x) 



ix^ 



1 - 



X In 



{u+ - u_) (l 

— M_ + U^U^ — 1 

— lt_ — U^U^ + 1 



(77) 



With y = l-^J f , and u± = 



1. Real part: collective enhancement factor 

The real part of the collectivity parameter governs the 
short-time MWA rate, and is given by 



dx \p{x)\ g{v,T,x). 



(78) 



Expanding to leading order in 1/t; gives the large- volume 
limit 

Mmf, = 3J^ Tdx^g^. (79) 

i;^oo V W Jo VJ^^ + 4x2 

We see that in this limit, the collective gain times the 
square-root of the volume, \/vfR, is a universal parame- 
ter, depending only on the Fresnel number, J- . Expand- 
ing to leading order in \/F gives the ^ 1 limiting 



lim fn — — ^_ , 

Vv:fJo 



dx X \p{x)\^ 



75 



(80) 



which gives the standard scaling oi D ^ nX^L for the 
optical depth, D = JrN. In the opposite limit of a 
very small Fresnel number, ^ 1, which corresponds 
necessarily to a large aspect ratio, fn has the limiting 
form 



lim fR=l\ - r dx \p{x)\^ = (81) 
.^c^-'^ 2V w io 14V^ 

which gives the expected scaling oi D ^ nX^L^, where 
Ld = W^/X is the diffraction length. 

The two limits can be merged via the Fade approxi- 
mant 



Ir ~ Co 



1 + CiJ" 



V 1 + C2F + C3 ' 



(82) 



with Co = 3.3660, ci = 0.1256, C2 = 0.3390, and cg = 
4.510 X 10^^, determined from the asymptotic formu- 
las and numerical fitting in the cross-over regime, which 
gives excellent agreement with the exact expression ( 78 1 . 
This is illustrated in Fig. [sj where we plot i/v/^ versus 
F for V = 00. Figure |3^ shows a wide range of J^- values 
on a log-log scale, whereas Fig. [3|d shows the crossover 
region on a linear scale. We note that y/vfji is effectively 
independent of v ioi v ^ 10^, which is generally satisfied 
for typical condensates containing N ^ 10® atoms. 

We see that for fixed condensate volume, the coeffi- 
cient fn, and therefore the optical depth for collective 
emission, D — f^N, is maximized at J- = 3.181, cor- 
responding to an aspect ratio of .4 = 0.420u^/''. For 
V ~ 10^, this gives A = 23.6 as the optimal aspect ra- 
tio with respect to the collectivity of the MWA process. 
For N = 10®, the corresponding maximum optical depth 
is then Dmax = 1050. For an aspect ratio of .4 = 10, 
at u = 10^, we have instead T = 10, and D = 853. 
For spherical geometry, we have instead J- — 215 and 
D = 202. 



2. Imaginary part: cross-phase modulation 

The imaginary part of the collectivity parameter is pro- 
portional to the induced dipole-dipole interaction energy 
of the sample. From a nonlinear optics perspective, we 
see from Eqs. (62), (63 1, and (64), this energy manifests 



itself as a source of cross-phase modulation between the 
ground and excited states. The magnitude of the imagi- 
nary part of / is given by 



= — dx\Pixfh{v,:F,x). 



(83) 
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FIG. 3: (Color online) The real part of cooperative parameter, 
fn, as a function of the Fresnel number, T. Discrete points 
are numerical results for the Thomas- Fermi profile ( 72 I for 



finite volumes, while the solid line corresponds to numerical 
results for v = <x. The (red) circles, (black) squares, and 
(blue) triangles correspond to v = 10^, v = 10*, and v = 10' , 
respectively. 



FIG. 4: (Golor online) The imaginary part of cooperative 
parameter, //, as a function of the Fresnel number, J^. The 
discrete points are numerical results for the Thomas-Fermi 
profile ( |72[ ) for finite volumes, while the solid line corresponds 
to the numerical result for v — oo. The (red)circles, (black) 
squares, and (blue) triangles correspond to v — 10^, v = 10**, 
and V = 10^, respectively. 



Expanding the integrand to leading order in povifers of 
1/v gives, 



lim // 



6 T 



dxx \p{x)\ 



In 



I" 



(84) 

Again, we see that y/vfj is a universal parameter, de- 
pending only on !F in the large- volume limit. Taking the 
limit — > CXI gives 



12 



lim fi ~ 



dx x^ \p{x)\'' 



180 1 



(85) 

scaling therefore as LlW"^. Similarly, taking the limit 
— >■ 0, corresponding to a very large aspect ratio, gives 



dx \p{xy\' In — 



lim // = 

15J" /1517 



Jy/^ V 1260 



-7e + ln2-lnJ" ,(86) 



whose pre-factor scales as 1/L, but which diverges loga- 
rithmically with decreasing 

The two limits can be joined by the Pade approximant. 



(^0-1"!^) (1 + ^1-^) 
7^AJ l + d2T + d3T'^ + dAT^' 



(87) 



where do = 1-453, di = 0.09402, da = 1.223, d^ = 0.5531, 
d4 = 0.1215, and = 0.01138. This is vahdated in Fig. 
[4] where we plot y/vfi versus for ?; = 10'', 10^, 10^. 
Figure |4^ shows a wide range of J^- values on a log-log 
scale, while Fig. |4Jd shows the cross-over region on a 
linear scale. 

The strength of cross-phase modulation relative to the 
collective gain is given hy x = fi / fn- In the limit <C 
1, corresponding to a highly-elongated condensate, the 
asymptotic form of this parameter is 

X w 27r(1.453 -ln J") > 1. 

In the opposite regime, J-" 3> 1, we have instead 

24 



(89) 
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FIG. 5: (Color online) The nonlinear coefficient x ~ fil Ir 
versus Fresnel number, on a log-log scale. The solid line 
gives the numerical result for w = oo, and the red dot indicates 
the location of the maximum /_r at = 3.18. 



which is small provided ^ 8. The full behavior of % 
is shown in Fi g. [5} where we plot the ratio of the Pade 
approximants (87) and (82 1 versus F. We see that x 



decreases monotonically with increasing which corre- 
sponds to a decreasing aspect ratio. 



C. QCM dynamics: analytic approximations 

In this section, we study the dynamics of the Raman 



MWA system by solving the QCM equations (62)-(|64[), 



neglecting atom losses from the pump and signal mode 
space due to spontaneous emission. In the absence of 
losses, the QCM equations become 



dtCl = Z^Ce + ^ (1 + ZX) |Ce|^Ci 

dtCs = ^ (1 + ix) IcepCs 



(90) 
(91) 



" J(1-*X) E Ic.fce.(92) 



We will now derive approximate analytic solutions for 
both the far-detuned Stark and near-resonance Rabi- 
regime. 



1. Stark regime 

In the Stark regime, defined as \ fl\'^+r^D^ < the 
excited state amplitude adiabatically follows the pump- 
matter-wave amplitude. For negligible losses, we than 

,P w 1. The adiabatic solu- 



have |c2p « and |cip 
tion for Ce is then 



2A-X 



(93) 



which leads to 

drCi 



R 



iA 



:(l + *x)|c. 



ci, 



drCs = 



R , 



where 



R 



-ix) |ci| 



{2A-XY + 1 



(94) 
(95) 

(96) 



is the dimensionless MWA rate constant. Using |ci| 
1 — |csp, this leads to the rate equation 



drUs = -R(l - ns)ns 



(97) 



p is the fractional population of the signal 



where n< 
mode. 

This has the analytic solution. 



ns{t) 



71^(0)6 



R.t 



l + ns[Q) (e«*-l)' 
where R is the MWA rate, given in original units by 

K ' 



(98) 



(2A - xr£')2 + r2i:)2 ■ 



(99) 



This shows that the true resonance condition is A = 
Xr_D/2, corresponding to lol = We — Wi — \TD. This 
represents a density-dependent red-shift in the resonance 
due to the atomic dipole-dipole interaction. The FWHM 
of the resonance is FD, which is larger than the natural 
linewidth, F, due to superradiance-broadening. 

Expanding to leading-order in the detuning, we obtain 
the Stark-regime MWA rate. 



Rs 



\VLYTD 
4A2 ' 



(100) 



which is larger than the spontaneous photon-scattering 
rate by a factor of the optical depth, D, characteristic of 
a collective emission process. 

From Eq. (|98| we see that for small initial signal frac- 
tion, 71^(0) ^ 1, there are two distinct time-scales. For 
times short-enough such that ns{0)e^*- <C 1, depletion of 
the pump condensate is negligible. In this regime, we 
have 



ns{t) « ns{Q)e 



Rt 



(101) 



so that R is clearly identified as the MWA exponential 
gain-rate. For longer times, such that ns(0)e^* ^ 1, the 
pump condensate is depleted, resulting in the saturated 
response limf^oo "■s(i) = 1- 

The signal-mode phase-shift, defined via Cg — 
/7^g«0s(t)^ is governed by 



Or 



(102) 
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which yields the analytic solution 



ln(ni(0) +n,(0)e^*)] . (103) 



This represents a nonlinear phase-shift that can be at- 
tributed to the inter-atomic dipole-dipole interaction. In 
the limit t — >■ oo, we find the net MWA phase shift is 
given by 



lim (l)s{t) = ^In^rrr. 



(104) 



The sign of the phase-shift is positive, with a magnitude 
proportional to the log of the initial signal-mode pop- 
ulation fraction, hence a smaller seed-value results in a 
larger phase-shift. This is readily interpreted as being 
due to a longer amplification time during which to ac- 
cumulate the phase shift. In the weak pumping regime, 
< (2A + + 1, the Hmiting value of the MWA 
phase-shift is independent of the detuning and pump in- 
tensity, depending only on the dimensionless parameter 
X- 



2. Rabi Regime 

For the case fi^ ^ (A — x/2)^, the atomic transition 
will be saturated, and the the atoms will exhibit Rabi- 
oscillations between states |1) and |e). To develop ana- 
lytic approximations for this regime, we first introduce 
the dressed-state amplitudes 



1 



c± = -^(ci±Ce)eT^"-/2, 



(105) 



which satisfy the initial conditions c±(0) = ■\/ni(0)/2. 
We can construct secular equations of motion for c± (t) 
by dropping all fast-rotating terms, yielding 



drC± 



A I I 2 

-2^c± - -\c,\ c± 



\c±?)c±, 



1 



(i+*x)(ic+p + ic_nc. 



Taking into account the conservation law |c+| 



1, this leads directly to 



1 



(1 - ns)ns, 



and 



dAs = |(1 - n^), 



(106) 
(107) 

(108) 
(109) 



Which reproduce the results of Eqs. (981, (103), and 



instead by 



(|104|), but with the Rabi-regime MWA rate being given 

YD , , 



This shows that the signal-mode amplitude and phase 
dynamics in the QCM model are universal as the detun- 
ing is swept across resonance, governed by a single rate 
constant. 



R = 



(2A-xri:')2-hr2D2-^2|r]2| 



(111) 



3. Numerics and Comparisons 



To verify the preceding an alysis, we numerically solve 
the mean-field equations (67 1- (68), for the case N = 10^ 
and V = 10^. For the Stark regime, we take A = 10, and 
(l = 1, and choose T = 10, corresponding to an aspect 
ratio of A = 10 and an optical depth oi D = 853. For 
the Rabi regime, we take A = 0, f2 = 10, and choose 

= 3.18, corresponding to an aspect ratio of A = 24 
and a maximized optical depth of -D = 1048. We consider 
two sets of initial conditions, ni{Q) — 0.9, ns{Q) — 0.1, 
ne(0) = 0, corresponding to a typical amplification ex- 
periment; and ni{0) — ng{0) = 0.5, n^{0) = 0, corre- 



sponding to a state-transfer experiment. 

In figure |6] we plot the signal population Ug (t) versus 
Rt, where R is given by Eq. ( 111 ). The solid lines show 
the universal analytic curve ( |98p for the two initial con- 
ditions, and the data points show numerical results cor- 
responding to Stark, Zeno, and Rabi regimes. In figure 
[7j we plot the signal-mode phase-shift acquired during 
MWA, comparing our analytic result to the exact numer- 
ical solution of Eqs. (67 1- (68 1. The figure plots (j>s{t)/x, 
which normalizes out the dependence on the Fresnel num- 
ber, J^, versus Rt, with R again given by ( 111 I. The solid 
lines show the universal analytic curve ( 103 1 for the two 



initial conditions, with the data points again showing nu- 
merical results corresponding to Stark, Zeno, and Rabi 
regimes. 

In conclusion, for a given MWA setup, the fastest am- 
plification of the signal wave is achieved in the Rabi 
regime, using an strong, resonant driving laser. Within 
the QCM approximation, both the population- and 
phase-dynamics is universal across the three main op- 
erating regimes. The population dynamics depends on 
geometry through the optical depth, D = fjiN, and is 
therefore optimized for speed by taking — 3.18. The 
MWA phase-shift approaches a limiting value that scales 
linearly with x = fi/fRj ^-^id is thus minimized by tak- 
ing as large as possible. In this study that corresponds 
to a spherical condensate, as we have not considered the 
case of a pancake BEC geometry. For v = 10^, going 
from J" = 3.17 (A = 23) to J" = 210 {A = 1.0) results 
in a decrease in the MWA gain rate i? by a factor of 5, 
with a decrease in the phase-shift by a factor 45. 
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Rt 



FIG. 6: (Color online) The signal population, ns{t), plotted 
against Rt, revealing the universality of the QCM population- 
dynamics. The upper solid line corresponds to Eq. ( |98[ l for 
the case ni(0) = ns{0) = 0.5, with the lower solid line showing 
the case ni(0) — 0.9, ns(0) = 0.1. The data points show 
the numerical results of the full QCM equations, (67l-(68l. 



The (blue) diamonds correspond to the Stark-regime, with 
A = lOFD, n = VD, and the (red) triangles show the Rabi- 
regime, with A = and = lOFD. 



f At) 




Rt 



FIG. 7: (Color online) The re-scaled signal-mode phase, 
4>s{t)/Xt plotted against Rt, revealing the universality of the 
QCM phase-dynamics. The upper solid line corresponds to 
Eq. (1031 for the case ni(0) = ns(0) = 0.5, with the lower 
solid line showing the case ni(0) = 0.9, ns(0) = 0.1. The data 
points show the numerical results of the full QCM equations, 
(67l-(68l. The (blue) diamonds correspond to the Stark- 



regime, with A = lOFD, f2 = TD, and the (red) triangles 
show the Rabi-regime, with A = and Q — lOFD. 



4- Mode-competition and the MWA threshold 

Atom losses in the QCM model are described by the 7 
term in Eq. (69 1. These losses correspond to the deple- 



tion of the atomic mean-field associated with the spon- 
taneous emission of incoherent radiation. In this section, 
we consider the possibility that some spontaneous pho- 
tons are emitted into modes that, while orthogonal to 
the idler mode, still lie within the solid-angle for end- 
fire emission, and therefore have a large gain-factor for 
amplification. Such amplified spontaneous emission is 



typically referred to in this situation as matter-wave su- 
perradiance [3]. Here we include atom losses in our an- 
alytical approximations, after scaling time to the MWA 
rate constant, R, the superradiant decay of the pump 
condensate is only weakly dependent on the number of 
end-fire modes. Furthermore, we show that the effect 
of pump depletion due to incoherent emission depends 
solely on the optical depth. 

The adiabatic elimination of the excited state in the 
Stark regime, and the introduction of dressed-states in 
the Rabi regime lead to two seemingly incongruous mean- 
field theories, given by Eqs. ( [MpS] ) and ( [l06p07| . With 



the introduction of the pump population, Np, given in 
the Stark regime by 



iVi and in the Rabi regime 
by Np = N+ + N_, i.e. the net population of the two 
dressed states defined by |105[ the population dynamics 
of the pump and signal modes obey a single set of rate 
equations. 



dtNp = -§NpNs 
dfN. = ?^.N„N, 



N 



(112) 
(113) 



These equations can be generalized to include incoherent 
emission and competing superradiance via the addition 
of spontaneous emission terms, and the inclusion of a 
complete set of atomic field modes, {|k, j)}, where k la- 
bels the spatial mode, and j the internal hyperfine level, 
which yields 



dtNp 



R 

'n 



Np[n, + 1 + J2 I (114) 



R 



dtN, = -NpiN, + 1), 
R 

dtN^j = ^iVpSkj(^kj +1), 



where 



fk] 

r fR 



(115) 
(116) 

(117) 



is the ratio ratio of the gain-factor for the |j, k) mode to 
that of the signal mode, satisfying 



1 



(118) 



For the level schemes depicted in figures [T] and [2j because 
the upper and lower levels are both J = 1 levels, decay 
to the mj — ground state is forbidden. Thus there 
are two decay channels with Fi = r2 = F, so that 7 = 
2F. Henceforth we shall assume these conditions hold. 
Introducing the fractional populations Up, n^, and ntj, 
and scaling time by the rate constant R, we arrive at 



dtUp 




(119) 
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1 

N ) ' 



N 



(120) 
(121) 



We note that the incoherent loss rate in (119) of 2/D, 



derived from these ad-hoc rate equations is exactly the 



rate derived from the full quantum-field theory ( 69 1 , 7 = 
2/D. 

For a typical elliptical condensate with 1 we have 



l-|d,f 



V cos^ Ok + sin^ 0]^ ' 



(122) 



where 0k is angle between k and z. The denominator 
gives the effect of geometry, whereas the numerator gives 
the dipole-emission pattern. For condensates with a large 
aspect ratio, >• 1, the transition to a strongly sup- 
pressed gain occurs at the end-fire angle Oef = W/L. 
Thus, all side modes can be classified into end-fire (EF) 
modes, corresponding to < (^ef, and non-end-fire 
(NEF) modes, corresponding to 0k > 0ef- We will as- 
sume that a given NEF mode never builds up any popu- 
lation, but that EF modes may become macroscopically 
populated due to superradiance. To greatly simplify the 
analytics, we can therefore approximate Skj by a step 
function 



Ski ~ H'^EF 



(123) 



Using the conservation law Hp + Us + X^kj ""kj = 1 to 
eliminate ns{t), results in a closed equation for rip, 



dfUp 



-Up [1 - Up + — 



which has the solution 



ip{t) = 



np{0) (1 + I,) 



np{0) 



1 



np(0)) e^*(i+2/-D) ■ 



(124) 



(125) 



In the limit D <C 1, this describes ordinary exponential 
decay, 



lim n„(t) = n„(0)e-2«/^ 



(126) 



while in the opposite limit of a large optical depth, it 
describes super-exponential decay 



lim np{t) = — 

D— ^00 1 



np{0)e 



-Rt 



!)■ 



(127) 



More important to consider is the signal mode popu- 
lation, given in the limit iV ^ 1 by 



ns(t) 



«s(0)(l + ^: 



(0) + -| + (l-n,(0))e-«*(i+2/^)' 



(128) 



To see the existence of an MWA threshold, consider the 
final population of the signal mode. 



,(00) 



.(0)(1 



D/ 



2_ 
D 



(129) 



In the limit Z) — >■ 0, this gives lim£)_>.o ns(oo) — ris(O), 
so that no MWA has occurred. In the opposite limit 
D — > CXI, it gives Hmu^oo 't-s(oo) — 1, which is the result 
with no competing incoherent losses. 

The threshold for MWA can therefore be defined as 
ns{oo) = 0.5, which means that lossless amplification of 
an initial signal of using a pump condensate of Np 
atoms requires an optical depth oi D ^ Dc, where 



Np 



(130) 



In our numerical simulations, A"p(0) = 9 x 10'', Ns{0) — 
1 X 10^, so that Dc = 16. Choosing = 3.18 to maximize 
the optical depth results in D — 1050, while a typical as- 
pect of = 10 results in D ~ 853, so that losses are 
essentially negligible during the MWA process. From a 
different perspective, we can conclude that lossless ampli- 
fication of a signal condensate of Ng atoms using a pump 
condensate of Np atoms having optical depth D requires 



2 + D 



(131) 



For Np ^ 10^ and D ~ 10'^, the initial signal popula- 
tion must therefore satisfy Ns » lO'^, i.e. rts(O) ^ 10~^, 
for strong amplification to occur. The final occupation 
of a typical competing EF mode is given for D ^ 1 by 
Nef{'^) — Np/Ns > 1. For the parameters of our simu- 
lations there are M ~ JF^ EF modes, so the total num- 
ber of atoms transferred to competing modes for the case 
J" = 3.18 is MNef ^ 10^ which is 0.01% of the total 
atoms. For comparison, about 2% of the total atoms are 
lost {nioss = 1 — ^sioo)), the vast majority of which are 
emitted into the NEF modes, which vastly outnumber 
the EF modes. 



IV. SPATIAL DEPENDENCE OF THE 
MATTER- WAVE DYNAMICS 

The primary limitation of the QCM model is not that 
it neglects the spatial dependence of the light fields inside 
the condensates, e.g. the idler field is zero at the front of 
the BEC, and strongest towards the back. The limitation 
is that it does not allow for a similar spatial dependence 
of thewavefunction of the condensate, which in the case 
of Rayleigh MWA would correspond to a spatial modula- 
tion in the visibility of the matter-wave density grating. 
In practice, such effects have been shown to account for 
several important experimental observations [El l40ti42] . 

To incorporate these spatial propagation effects, we 
now turn to numerical simulation of the mean-field Equa- 
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tions (36)- (38 1. Assuming a cylindrically symmetric sys- 



tem, we can make the change of variables 



(132) 
(133) 



V'3(r) = ViVVie(p,^)e^'"^ 



the mean-field equations become 



(134) 



9rV'l(P, Z) = 
drtpeiP, z) = 

with 

G{p, z; p' , z') 



2^*V'e(p,2;) + - / p' dp' dz' G* {p, z; p' , z' )%IjI{p' , z')'ijji{p' , z')il}e{p, z,t), 



p dp dz G*{p,z;p ,z)'tpl{p ,z )ipsip ,z )ipe{p,z), 



-iA- 



3 



(135) 
(136) 



Mp,z) + -n^i{p,z) + - J2 / p'dp'dz'Gip,z;p',z')4>*{p',z')^eip',z')i^j{p,z), (137) 



1 ll giK(0')-fcl-(2-z')] 



e(0') 



(138) 



r 



where 

^(0') = kL^/{z - z'Y +p2 + p/2 _ 2pp'cos<^', (139) 

and as with the QCM model we have r = FDi, — 
n/TD, A = A/ri), with D = JrN. 

A. Adiabatic elimination in the Stark regime 

For the Stark regime, A'^ ^ |f2p -I- 1, we can adiabati- 
cally eliminate the excited state dynamics by expanding 
tipe to second-order in powers of 1/A. This gives 



n 

4A2 



V'e « 7^ («7 + 2A - 
which leads to 



(140) 



2 4A2 



(2A 



i^+iG+G*)[rii'i])i^ 



+ '-^3^G[rji]^p, 

2 4A2 ^^^^-^J^^ 



(141) 

(142) 



The first term in Eq. ( 141 ) describes the decay of the 



mean-field due to spontaneous emission. The second 
term is the usual spatially homogeneous AC-Stark shift, 
with the third term describing the spatial- variation to the 
AC-Stark shift due to any modulation of the pump laser 
amplitude caused by the atom-field interaction. The re- 
maining terms in Eqs. ( 141 1 and ( 142 ) describe the trans- 



fer of atoms from the pump-mode to the signal mode 
which lies at the heart of the MWA process. 



With the change of variables, r 
'016*^'^ , these simplify to 



4A^ 



\t' and Tpi 



dr'ips = 



-|^i + ^(G + G*)K*^;K 

G[rsi''i]i's (143) 

G*K>,]v;;. (144) 



Note that these equations are invariant under the trans- 
formation ips — >■ e~*'^^s, which means that the popu- 
lation and phase-shift dynamics are independent of any 
spatially uniform initial relative phase between and 
"05 ■ This shows that in the Stark-regime, the signal mode 
population and phase dynamics is unaffected by the spa- 
tially homogeneous part of the AC-Stark shift, and its 
dependence on the pump intensity and detuning disap- 
pears upon scaling time to the rate constant R — ^^f^^ , 
just as with the QCM model. 



B. Numerical simulations 

In this section we give results of numerical simulations 
of the mean-field theory given by Eqs. ( 135 1-( 1371. We 



will compare results using the full Green's function of Eq. 
(1381, to those using the Paraxial Green's function, Eq. 
(511, and the Longitudinal Green's function, Eq. (54 1 as 



well as the results of the QCM model, Eqs. (67)-(69 1. To 



study the limit of very large detuning, we instead solve 
the adiabatic equations, (143) and (144 1. The spatial 



mean-field equations are solved using a grid of 20 points 
along the radial direction and 200 points along the lon- 
gitudinal axis. This requires first computing and storing 
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the greens function as a 4000 x 4000 array, which pre- 
vents us from increasing the number of mesh points due 
to computer memory hmitations. One test of the ade- 
quacy of this grid is to use it to compute numerically the 
collectivity parameters and //. We find that for the 
full Green's function of Eq. (138), the deviation from 



the exact values of ffj and // is 1%. Based on this, we 
conclude that the 20 x 200 grid is sufficient to accurately 
model the true electric field inside the condensate during 
the MWA process. 

We consider a condensate of = 10^ atoms con- 
fined to a dimensionless volume of w = 10''. For initial 
conditions, we take rpi{r,0) = y'ni(0)(/>(r), '(/'^(r, 0) = 
•\/n^(O)0(r), and V'e(r, = 0, with ns(0) being the ini- 
tial signal population fraction, and ni(0) being the initial 
pump fraction, so that ni(0) -1-77.^(0) = 1, and with the 
mode function (f){r) given by the Thomas-Fermi wave- 
function of Eq. (72). The pump laser is taken to be 
spatially uniform, with a sudden turn-on at t = 0. 



1. Stark Regime 



For the off-resonant Stark regime, defined as ^ 
j-j2 _|_ p2£)2^ have numerically solved the adiabatic 
equations ( |143 ) and ( 144 ) for a variety of geometries and 
initial conditions. In Figure |8] we the signal-mode popu- 
lation fraction, ng(t) versus Rt, where 



flsit) 



dr\tp{r±, z) 



(145) 



and R is given by Eq. (111). The initial conditions are 
taken as ns(0) — 0.1 and ni(0) = 0.9, corresponding to a 
small-signal amplification process. The solid (black) line 
gives the results from the QCM (single-mode) model, the 
dash-dotted line (blue) shows the results from the longi- 
tudinal Green's function, the dashed (green) line shows 
the results from the Paraxial Green's function, and the 
dotted (red) line shows the results from the exact Green's 
function. 

In Figure [8^, we show the results for T = 3.18, cor- 
responding to an aspect ratio oi A — 24, which gives 
the maximum gain within the QCM model. In this case, 
we see that there is excellent agreement between the full 
and Paraxial Green's functions, which is expected for the 
case of a large aspect ratio. We also see that there is ex- 
act agreement between the single-mode approximation 
(QCM) and the multi-mode exact model at short times 
At longer times, the exact model saturates with only 
about half of the pump atoms transferred to the signal 
mode, whereas the single-mode theory predicts a com- 
plete transfer. The results from the longitudinal Green's 
function do not agree with either the QCM or exact mod- 
els, greatly overestimating the initial gain and final sat- 
uration level. 

In Figure[8]3, we depict the case F = 10, corresponding 
to a shorter aspect ratio of ^ = 10, which is typically of 
many MWA and superradiance experiments. Here we see 



that for short times there is again agreement between the 
single-mode QCM model and the exact model. At longer 
times, the paraxial approximation is not quite as good as 
with the A — 24 case. The inaccuracy of the longitudinal 
model is slightly improved. The final saturation level 
of the signal-mode population is seen to increase with 
respect to the = 3.18 case. 

Lastly, in Figure [8];, we show results for = 30, corre- 
sponding to an even shorter aspect ratio of ^ = 4.4. In 
this case, we now see good agreement between the Lon- 
gitudinal and exact Green's functions. Again the agree- 
ment with the QCM model is exact at short times, fol- 
lowed by a break-away and transition to a very slow rate 
of increase in the signal-mode population. 

We now turn our attention to the phase-dynamics of 
the signal matter- wave. Ideally, a MWA device should 
leave the signal-mode phase unaltered. From the QCM 
model, we expect that this is not possible within the cur- 
rent paradigm. If the induced phase shift can be calcu- 
lated precisely, however, then it can be taken into account 
so that an MWA device can still be employed for interfer- 
ometric applications. In figure[9] we plot the signal-mode 
phase divided by x versus Rt for the same parameters as 
in Figure |8] The signal mode phase has been spatially 
averaged, weighted by the signal-mode density, according 
to 



ns{t) 



dr|Vs(r,t)|2arg[0,(i)]. 



(146) 



From Fig. 8, we see that with respect to the phase dy- 
namics, the agreement between the QCM model and the 
exact result, while still perfect at very short times, di- 
verges much more quickly than the population dynamics, 
more so with increasing J-. Thus the closest agreement 
with the QCM model is in the T = 3.18 case. Note 
that the Longitudinal Green's function model predicts 
zero phase-shift in all cases. Thus even for the case of 
a large Fresnel number, where the Longitudinal model 
accurately predicts the population dynamics, as seen in 
Fig. Isfc, it does not accurately model the phase dynamics. 



2. Rabi Regime 

For the on-resonance Rabi regime, w e tak e A = = and 
Q, — 10, and solve numerically Eqs. (135)-(137). The 



signal-mode population dynamics is poltted versus Rt 
for initial conditions ns(0) = 0.1 and rii(O) = 0.9, in Fig. 
[To] FigureifTO^, b, and c correspond to = 3.18, 10, 
and 30, respectively. As with the Stark regime, we see 
good agreement between the full Green's function and the 
Paraxial Green's function for the highly-elongated F — 
3.18 case, with large discrepancies occurring for larger 
Fresnel numbers. In the opposite limit of a relatively 
small aspect ratio, we see for = 30 that there is good 
agreement instead between the full Green's function and 
the Longitudinal approximation. 
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FIG. 8: (Color online) The signal mode population frac- 
tion ns{t) versus Rt, for the Stark regime, determined by 
Eqs. (1431 and (1441, for the initial conditions Us — 0.1 and 
ni = 0.9, with iV = 10*^ and t; = lO'^. In all figures, the (black) 
solid line shows the results of the QCM model, the (red) dot- 
ted line shows the results for the full Green's function, the 
(green) dashed line shows the results for the paraxial Green's 
function, and the (blue) dash-dotted line shows the results 
for the Longitudinal Green's function. Figure (a) shows the 
case T = 3.18, corresponding to an aspect ratio of ^ = 24, 
figure (b) shows the case J- = 10.0, corresponding to an as- 
pect ratio of ^ = 10, and figure (c) shows the case — 30, 
corresponding to ^ = 4.4. 



In figure [TOp ., we see a very important contrast be- 
tween the on-resonant and the ofF-resonant dynamics in 
the maximum gain, T = 3.18 case. Unlike the off- 
resonant case of Fig. [8^, where the signal mode popula- 
tion determined via the full Green's function saturates at 
ns{t) ~ .5, well below the prediction of the longitudinal 
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FIG. 9: (Color online) The spatially averaged signal-mode 
phase, <^s(t) scaled by x s-nd plotted versus Rt for the Stark 
regime. All parameters are as in Figure |8] 



Green's function, in the on-resonant case, the saturated 
value of ng{l) is even slightly above that of the Longitu- 
dinal model. This effect is also seen in figures [TOjj and c, 
albeit the differences are less pronounced. 



3. Conclusions 

In figure [8j we see a striking difference between the 
prediction of the Longitudinal Green's function and both 
the exact and Paraxial Green's function, with the MWA 
process being strongly suppressed in the latter case. In 
Fig. [TOj we that on-resonance this suppression goes away. 
The difference between the off-resonant and on-resonant 
population dynamics can be attributed to the spatial- 
variation in the AC Stark-shift of the pump matter-wave 
which arises due to laser-depletion. The right-circularly 
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FIG. 10: (Color online) The signal mode population fraction, 
na{t) versus Rt for the Rabi regime, with A = and Q,l = 
10r_D. The initial conditions and geometric parameters are 
the same as in figure [8] 



polarized light field can be decomposed into a superpo- 
sition of the bare (undepleted) laser field, and a "de- 
pletion field", which is opposite in sign, as 57+(r,t) = 
+ ^^d(r, i). In the case of a large detuning, the AC 
Stark-shift of state 1 1) is given by 



^Stark 



l»+(r,^)P 

4A 
|»Lp 1 
4A 4A 
l^d(r,t)p 
4A 



{niQdir, t) + c.c) 



(147) 
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FIG. 11: (Color online) The spatially-averaged signal-mode 
phase scaled by x plotted versus Rt for the Rabi regime, 
with A = and Ql = lOFD.The initial conditions and geo- 
metric parameters are the same as in figure [8] 



The depletion field is generated by the induced polar- 
ization, and therefore scales as ftd ^ f^i/A. Thus to 
second-order in 1/A, we have 



A 



Stark 



4A 



8A2 



(148) 

This shows that the spatially-dependent part of the AC 
Stark-shift scales with the laser-intensity and detuning 
exactly as the MWA gain term. Thus the effect of pump- 
depletion is universal, i.e., it does not improve by going 
further off resonance and increasing the laser intensity. 
This means that for MWA, there is no such thing as the 
undepleted- pump approximation in the Stark regime. The 
MWA suppression due to the depletion field should be 
largest when the real-part of G is strongest, corresponds 
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FIG. 12: (Color online) Numerical results for the same pa- 
rameters as figure [sja), but with the spatially-dependant AC 
Stark term removed from the field equations, verifying that 
this removes the MWA supression. 



to the largest x^ ^^'^ hence the smallest as can be seen 
in the numerical simulations. 

To verify that this explanation is correct, we repro- 
duce Fig. [8^, but with the spatially dependent Stark- 
shift term removed, i.e. the second-term on the r.h.s. of 

where we 



Eq. (143). The result is shown in figure 12 



see that the effect disappears. To better understand the 
mechanism for MWA suppression in the Starkj regime, 
note that If the phase differential across the signal-mode 
is of order tt, then the idler field emitted at the front of 
the condensate will be reabsorbed rather than amplified. 
That this is indeed what is going on is shown in Fig. |13[ 
where we plot the field signal strength |'0sP and phase, 
aigips, at the center of the condensate, p = 0, as a func- 
tion of z and t. We see that initially the signal mode 
is built up towards the back of the BEC, as expected, 
but as the amplification wave moves towards the front, 
the signal initially built up at the back is reabsorbed. 
Thus the signal amplification occurs not at the rear, as 
has been observed experimentally, but at the front of the 
BEC. The reason this effect was not seen experimentally, 
is that previous experiments pump from the side, so that 
laser depletion, and the associated AC Stark modulation 
does not play a significant role. 

Even in the absence of this re-absorption effect, the 
fully spatial models show a cross-over from a fast-gain 
process, accurately captured by the QCM model, to 
a slow-gain process. This cross-over is present in all 
regimes, and is not related to the spatial AC Stark-shift. 
Instead, it is related to the spatial build-up of the idler 
field, so that atoms at the front of the condensate see 
a very weak field and are inefficiently transferred from 
the pump to signal modes, as described by previous au- 
thors. Note that even with this effect included, there 
is still exact agreement with the QCM model at short 
times. This is because the spatial structure of the idler 
and pump fields is exactly included in the calculation of 
fa and //, the parameters that govern the QCM model. 
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FIG. 13: (Color online) Signal mode density at p = ver- 
sus Rt and kLZ, for the Stark regime with — 3.18 and 
ns(0) = 0.1. The shading indicates the phase of the signal 
amplitude. In (a) we see that the signal mode is initially 
amplified at the back of the condensate, but is then reab- 
sorbed, so that signal-mode amplification occurs at the front 
of the BEC. Note that this effect would not occur when the 
pump laser is perpendicular to the long-axis. In (b) we show 
the numerical results with the spatially dependent AC-Stark 
term removed, showing that reabsorption does not occur, so 
that the amplification is concentrated towards the back of the 
condensate. 



It is only the spatially-dependent response of the matter- 
waves that is not included. The cross-over into this slow 
gain regime is due to the fact that the signal matter- 
wave is concentrated mainly in the back of the initial 
condensate region, and the source mode is concentrated 
in the front region due to depletion at the back. These 
two modes become highly mis-matched, so that the idler 
field generated by coherence between them becomes very 
weak. 

The single-mode QCM model predicted that the dy- 
namics in the Stark and Rabi regimes is universal, with 
the only difference being in the underlying time-scale of 
the dynamics. In the full multi-mode model, we see that 
in addition to a faster dynamics, the on-resonant Rabi 
regime has significantly better performance due to the 
absence of the spatially dependent AC Stark-shift of the 
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laser-depletion field. For all regimes, we see that the 
signal- mode phase-shift scales as Xi which is purely a 
function of the Fresnel number T . By decreasing the 
Fresnel number, one can increase the MWA transfer rate, 
at the expense of increasing also the non-linear phase- 
shift. Conversely, by increasing the Fresnel number, one 
can reduce the phase-shift at the cost of reducing the 
MWA transfer rate. For a given choice of T ^ and hence 
X, it is always advantageous in terms of time to work in 
the on-resonant Rabi regime. Furthermore, if one choses 
to maximize the MWA gain by taking T « 3.18, then 
the Rabi-regime has the additional advantage of the ab- 
sence of idler reabsorption due to the AC Stark-shift of 
the laser depletion field. 



C. Gain versus losses 

Our numerical simulations include atom losses, i.e. de- 
cay of the atomic mean-field intensity, which occur due 
to incoherent scattering of the laser photons into modes 
other than the idler mode. In this process, a recoiling 
atom is entangled with a scattered photon, so that both 
the atomic and optical mean-fields are decreased. This 
process is suppressed relative to idler emission by a fac- 
tor of the optical depth D = /^iV, due to the collective 
enhancement of the latter process. Thus losses are mini- 
mized by maximizing Z?, which occurs at = 3.18. 

One characterization of the efficiency of the MWA pro- 
cess is the gain-to-loss ratio. In Fig. ( 14 ) we plot the gain 



0.015 



versus loss trajectory of the signal mode for both the 
on-resonant Rabi regime (thick red lines), and the off- 
resonant Stark regime (thin blue lines). The atom loss 
percentage is defined as Loss = ns{t) + ni{t) +rte(<), and 
the gain factor is defined by Gain = {ns{t)—ns{0))/ns{0). 
The dashed lines correspond to = 3.18, the dot-dashed 
lines to = 10, and the dotted lines to J" = 30. We see 
that by this metric, resonant MWA with = 3.18 is 
optimal. It is fortuitous that this also corresponds to 
the case where the MWA transfer process occurs on the 
fastest time-scale. 



D. State transfer dynamics 

Another potential application of the MWA system is 
the transfer of a quantum state from a matter-wave field 
onto an optical field. For example, the relative phase 
between the left and right modes of a double-well BEC 
could be measured by mapping the relative phase onto an 
idler field, which is then measured via homodyne detec- 
tion with the pump laser. This would be accomplished 
by first removing the double-well potential barrier, so 
that the two condensates move towards and through each 
other, and at the moment in which the two modes over- 
lap, a single pump laser is turned on to amplify on of 
the modes at the expense of the other. The phase of the 
idler field generated by this process will be given by the 
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FIG. 14: (Color online) Here we plot the atom-loss percentage 
versus the signal-mode gain factor calcuated using the full 
Green's function. The (red) thick lines correspond to the 
on-resonant Rabi regime, and the (blue) thin lines to the ofF- 
resonant Stark regime. The dashed lines correspond to 7^ = 
3.18, the dash-dotted lines to = 10.0, and the dotted lines 
to = 30.0, we see that the lowest loss-to-gain ratio occurs 
on-resonance, and under conditions of maximum gain, = 
3.18. 



relative phase between the two condensates. In mean- 
field theory, we cannot make claims regarding squeezed 
states, however, we believe that conservation laws will be 
sufhcient to guarantee that number- and phase-squeezed 
states of the two-mode condensate can be mapped onto 
the idler pulse without destroying the squeezing. 

To simulate this process, we again perform numeri- 
cal simulations, but now with initial conditions ^^(O) = 
ni(0) = 0.5. As was seen in the previous section, we find 
that complete transfer does not occur, due to the fact 
that atom transfer is inefficient towards the front of the 
BEC, where the idler field is very weak. Even for a 50% 
transfer, the number of photons scattered into the idler 
beam will be of order iV/4, so that the shot-noise of the 
idler pulse will be very close to that of the initial par of 
condensates. In Figures [T5| and [T6| we plot the idler field 
at the center of the condensate, r_L = 0, versus Rt along 
the X-axis and k^z along the y-axis. The height of the 
surface corresponds to the idler field intensity and the 
color corresponds to its phase. 

In Figure[T5]we show the results for the far off-resonant 
Stark regime, where the idler field is defined as 



(149) 



Figures (a), (b), and (c) correspond to the Fresnel num- 
bers J- — 3.18, 10, and 30, respectively. For T = 3.18, 
we see a small variation in the idler phase over space 
and time, which can be attributed to the spatial varia- 
tion in the AC-Stark shift of the source matter-wave, as 
described in Sec. IIVB 31 That this effect becomes less 
visible as is increased is due to the fact that this phase- 
shift is proportional to Xj which decreases with increas- 
ing J^. The second interesting thing that we see in the 
plot is that at long times, there is a relatively weak idler 
field inside the condensate that does not propagate out of 
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FIG. 15: (Color online) Idler field intensity at p = versus Rt 
and k^z for the off-resonant Stark regime, for the initial condi- 
tions ns{0) = ni(0) = 0.5, corresponding to a matter-to-light 
state transfer protocol. Figures (a), (b), and (c) correspond 
to T = 3.18, 10.0, and 30.0, respectively. 



FIG. 16: (Color online) Idler field intensity at p = versus 
Rt and kLZ for the on-resonant Rabi regime, with A = and 
JIl ~ lOF-D, and for the initial conditions ns{0) = ni(0) = 
0.5, corresponding to a matter-to-light state transfer protocol. 
Figures (a), (b), and (c) correspond to = 3.18, 10.0, and 
30.0, respectively. 



the condensate, but is instead re-absorbed at roughly the 
mid-point, z = 0. This is the idler field during the sec- 
ond, much slower stage of the MWA process, where the 
signal and source condensates have evolved into spatially 
mismatched modes. Note that the region at large t and 
z where the idler phase increases towards tt corresponds 
to a region where the idler amplitude is very small. 



In Figure 16 we show the results for the on-resonant 
Rabi regime, where the idler field is defined as 



ni{rj_, z,t) = G[-4)*^e]{r±,z,t). 



(150) 



Figures (a), (b), and (c) correspond to the Fresnel num- 
bers T = 3.18, 10, and 30, respectively. Here the phase 
is taken modulo tt, so that the discontinuous phase slips 
that occur during each Rabi oscillation are suppressed. 
We see that the signature of the Rabi oscillations is a 
modulation of the idler amplitude at the Rabi-frequency 
of the driving laser. Aside from this modulation, we 
see that the spatial and temporal uniformity of the idler 
phase increases as the Fresnel number increases, as with 
the off-resonant case. 
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V. CONCLUSIONS 

Starting from the standard Markovian reservoir the- 
ory, wc have derived a Green's function approach to the 
coherent mixing of matter- and light- waves, suitable for 
modelling the propagation of electric fields inside Bose- 
Einstein condensates. We have found that at the mean- 
field level, this approach reproduces exactly the field of an 
ensemble of ideal radiating dipoles, and is exactly equiva- 
lent to solving Maxwell's wave equation. We have studied 
two approximations to the full Green's function, (1) the 
paraxial wave equation, that includes diffraction effects, 
and (2) the longitudinal wave equation, that neglects the 
transverse component of the Laplacian. Whereas (2) has 
been used by previous aiithors, we find that it is inad- 
equate for the task. It does not include the nonlinear 
phase-shifts associated with the dipole-dipole interaction, 
which manifest themselves as the real-part of the Green's 
function. It also does not accurately model the superra- 
diant gain when the length of the condensate length be- 
comes comparable to the diffraction length. The paraxial 
approximation, on the other hand, can give accurate re- 
sults provided the aspect ratio is much larger than unity. 

To demonstrate the utility of this approach, we have 
extended the studies of matter-wave amplification into 
the resonant driving regime. Wc have considered a co- 
propagating Raman MWA scheme and compared the off- 
resonant Stark regime to the on-resonant Rabi regime. 
We find that in the off-resonance regime, MWA is sup- 
pressed by a spatially dependent AC Stark shift associ- 
ated with the laser depletion field. As a consequence, 
MWA is most efficient in the on resonance regime. We 
have verified that the rate of mean-field atom losses due 
to spontaneous emission (heating) scales with detuning 



and laser intensity in the same way as the MWA gain, 
so that the gain-to-loss ratio is the same both on and off 
resonance, and determined solely by the optical depth of 
the condensate. 

In addition, we have compared the single-mode approx- 
imation for the matter-waves to the fully spatial propa- 
gation, and find that the single-mode model gives ex- 
act agreement for the short-time evolution. The MWA 
gain is determed by the optical depth of the conden- 
sate. The optical depth D, defined in the short time 
regime, can therefore be computed analytically. We have 
found that upon scaling by y^, where v = k^V, with 
V being the volume, that for sufficiently large v the 
gain exhibits a universal behavior which depends only 
on the Frcsnel number T = k^W'^/L, which is the ra- 
tio of the diffraction length, Ld = kLW"^, to the BEG 
length, L. For L <C L^, we obtain the usual formula 
D ^ n}?L, but for the case i » L^, it saturates instead 
at f ~ nX'^Ld- This leads to a maximum optical depth 
of Djnax ~ N/Vv, occuring at L ~ L^, or more precisely, 
F = 3.18. In addition to the gain, there is also an as- 
sociated MWA phase shift, which scales as x, the ratio 
of the imaginary collectivity parameter to the real-part. 
Thus one can generalize the optical depth to a complex 
parameter Dcompiex = D{1 - ix)- 

To summarize, we propose that it is optimal to use 
resonant driving, instead of the commonly used off- 
resonance driving, to perform Raman MWA for the pur- 
poses of amplifying a matter wave and/or mapping quan- 
tum states between optical and matter waves. The rea- 
son being that the gain-to-noise ratio is slightly higher 
and the AG-Stark shift vanishes. This requires a short 
(~ Ins) pulse with a Rabi frequency satisfying flL > 
TD ~ 10^°. 
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